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ABSTRACT 

We have measured 492 redshifts (311 new) in the direction of the poor cluster 
AWM 7 and have identified 179 cluster members (73 new). We use two independent 
methods to derive a self-consistent mass profile, under the assumptions that the 
absorption-line galaxies are virialized and that they trace an underlying Navarro, 
Prenk & White (1997) dark matter profile: (1) we fit such an NFW profile to the radial 
distribution of galaxy positions and to the velocity dispersion profile; (2) we apply the 
virial mass estimator to the cluster. With these assumptions, the two independent 
mass estimates agree to ~15% within 1.7 h~'^ Mpc, the radial extent of our data; 
we find an enclosed mass ~ (3 ± 0.5) X 10^^ h-^ Mq. The largest potential source of 
systematic error is the inclusion of young emission-line galaxies in the mass estimate. 

We investigate the behavior of the surface term correction to the virial mass 
estimator under several assumptions about the velocity anisotropy profile, still within 
the context of the NFW model, and remark on the sensitivity of derived mass profiles 
to outliers. We find that one must have data out to a large radius in order to determine 
the mass robustly, and that the surface term correction is unreliable at small radii. 

Subject headings: galaxies: clusters: individual (AWM 7) — galaxies: kinematics and 
dynamics 



1. Introduction 



The properties of galaxy clusters play an important role in constraining cosmological 
parameters and models of large-scale structure formation. Because clusters are the most massive 
collapsed aggregates in the universe, the cluster mass function and its evolution with redshift 
constrains the spectrum of density fiuctuations (e.g. Eke, Cole & Frenk 1996, Bahcall &: Cen 
1992). Accurate cluster masses can also constrain the ratio of baryonic to total mass, and thus VIq 
(White et al. 1993). 



^ Observations reported in this paper were obtained at the Muhiple Mirror Telescope Observatory, a facihty 
operated jointly by the University of Arizona and the Smithsonian Institution; at the Whipple Observatory, a 
facility operated jointly by the Smithsonian Astrophysical Observatory and Harvard University, and at the WIYN 
Observatory, a joint facility of the University of Wisconsin-Madison, Indiana University, Yale University, and the 
National Optical Astronomy Observatories. 
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There are several standard ways of computing cluster masses. Optical redshift surveys yield a 

mass estimate through the application of the virial theorem, on the assumption that the galaxies 
trace the overall cluster mass distribution; this method was first applied to Coma by Zwicky 
(1933). More sophisticated variants of this method include a surface term to account for the 
unsampled portion of the cluster (The & White 1986), and use the Jeans equation to account for 
the effects of orbital anisotropy (e.g. Binney &; Tremaine 1987). For clusters with many hundreds 
of redshifts, the distribution of infalling galaxies in redshift space measures the gravitational 
potential of the dark matter halo (Diaferio & Gcller 1997); this method has only been successfully 
applied to Coma so far (Gcller, Diaferio & Kurtz 1999). Observations of hot intraclustcr gas in the 
X-ray provide information about the gravitational potential; the added assumption of hydrostatic 
equilibrium yields a mass (Bahcall & Sarazin 1977). Gravitational lensing of background galaxies 
can also be used to estimate cluster masses (Webster 1985, Bartelmann et al. 1996 and references 
therein). Comparing cluster masses derived by these independent methods can provide insight 
into the internal structure of clusters (Loeb &; Mao 1994) and provides a consistency check on the 
mass derived by any one method. 

While the distribution of total cluster masses constrains global cosmological parameters, 
knowledge of the density (or equivalently, mass) profiles of individual clusters constrains models 
of cluster formation. iV-body simulations of hierarchical clustering by Navarro, Frcnk & White 
(1995,1996,1997, hereafter NFW) suggest that a universal two-parameter density model can 
adequately describe clusters over a broad range of masses. Until recently, few clusters had enough 
measured redshifts to constrain the mass profile: the advent of multi-fiber spectrographs has now 
made it possible to measure several hundred redshifts in a cluster within a reasonable time. 

Our goal is to measure the mass profile of the cluster AWM 7, selected by Albert et al. 
(1977) on the basis of the cD galaxy at its center. Our sample consists of 179 cluster galaxies with 
measured redshifts; this work expands the sample of Koranyi et al. (1998) by 70%, extending the 
surveyed region to a projected radius of 1.7 Mpc, and reaching fainter magnitudes in the 
core. We confirm the cold core seen in Koranyi et al. (1998), and fit a NFW mass profile from 
the projected distribution of radii for comparison with the virial mass profile; we consider also 
the effect of orbital anisotropics (still assuming the NFW profile) on the mass estimates. The two 
profiles agree remarkably well within our restricted context; the assumption of the NFW model 
produces self-consistent mass estimates using these two methods. Relaxing the assumption that 
the galaxies trace an underlying NFW dark matter distribution would result in a larger range of 
admissible masses (cf. The Sz White 1986). 

In §2 and §3, we describe the observations and data reduction. In §4 we describe the data, and 
in §5 we justify the choice of the NFW model, compute the resultant mass profiles, and comment 
on the reliability of the surface term correction. We defer discussion of some details until §6, and 
conclude in §7. We use Hq = 100 h km s"^ Mpc~^ throughout. 
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2. Observations 

We observed a total of 593 distinct targets in the direction of AWM 7 from 1995 to 1999. We 
measured redshifts for 492 objects; 179 are cluster galaxies (with redshifts in the range 2500-7500 
km s~^), making AWM 7 one of the best-sampled clusters in the sky. Our sample extends over 
1.5 degrees in radius from the cluster center, and reaches magnitudes as faint as i? = 18 in the 
core. We collected data from three different telescopes, with different sky coverage and different 
magnitude limits; we describe the details of each in turn. In each case we chose the spectral 
coverage so that both the Ca H-K break at AA3933, 3968 A and Ha A6562 A would faU on the 
detector at the cluster redshift. These spectral features are important for determining redshifts by 
cross-correlation with template spectra. We list the redshifts for all measured cluster galaxies in 
Table |. 

2.1. FAST Observations 

We used the FAST Spectrograph (Fabricant et al. 1998) mounted on the 1.5-m Tillinghast 
telescope at the F.L. Whipple Observatory on Mt. Hopkins to observe 302 targets over ~10 
square degrees of sky centered on NGC 1129. We selected the targets from the Digitized Palomar 
Observatory Sky Survey scans. We used a 3 arcsecond slit and a 300 line mm~^ grating, resulting 
in a resolution of -^6 A and spectral coverage of 3600-7600 A. FAST can reliably measure redshifts 
for galaxies down to R ~ 15.5, and can measure redshifts for galaxies as faint as i? ~ 16 if the 
central surface brightness is sufficiently high. 159 FAST targets are new to this paper; 51 of these 
are cluster members. 

2.2. MMT Observations 

We observed 45 targets with the Blue Channel spectrograph of the MMT on the nights of 
1996 Dec. 5 and 6, using a 2 arcsecond slit and a 300 line mm~^ grating, resulting in a resolution 
of ~10.6 A and spectral coverage of 3000-8000 A. These targets are in the central 75' x 75' 
(1.15 X l.lbh^^ Mpc) of the cluster, and are complete to -R = 16.5 in that region. These data have 
been previously published in Koranyi et al. (1998). 

2.3. Hydra Observations 

We observed 250 targets (under the queue observing scheme) with the Hydra multi-fiber 
bench spectrograph mounted on the WIYN telescope at Kitt Peak on the nights of 1998 Dec. 20 
and 21; these data are new to this paper. There were 4 separate fiber configurations; 60 objects 
were observed in more than one configuration. There were either four or five 25-minute exposures 
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for each configuration. We used 3 arcsecond fibers, and a 400 line mm ^ grating, resulting in a 

resolution of ~7 A and spectral coverage of 3800-7000 A. The four configurations were observed 
under strongly varying conditions, resulting in dramatic diff'erences in throughput. Two of the 
configurations suffered from wind shake, and a third from cirrus. The remaining configuration 
was observed under stable, clear skies with good seeing and resulted in a redshift for almost every 
target. 

In total, we were able to measure redshifts for 154 of the 250 targets (152 new), sometimes 
by combining spectra of the same object observed under different configurations. Only 22 of these 
are cluster members. The Hydra targets are all within the central 1.4° x 1.4° of the cluster. 

The magnitude completeness of the Hydra sample is complicated. Not only were we unable 
to measure a redshift for every target, but the original target list itself did not include every 
galaxy (up to some limit) in a magnitude-ordered list because of physical constraints in allocating 
fibers to objects. Moreover, the initial magnitude ranking was not by total galaxy luminosity, but 
rather by luminosity within the central 3 arcseconds subtended by the fibers. We will discuss the 
photometry and luminosity function of AWM 7 in a later paper; here, we are interested primarily 
in the kinematics of the cluster and the incompleteness is not a limitation on the analysis. 

3. Data Reduction 

We reduced all spectra using the xcsao task of the RVSAO package in IRAF (Kurtz & 
Mink 1998). xcsao measures redshifts by cross-correlating the observed spectra against a suite 
of template spectra. The quality of the match is measured by the parameter R, an estimate of 
the ratio of the highest peak to the average noise peak in the correlation function. The template 
spectra are themselves empirically derived from an archive of FAST observations, and are thus 
strictly speaking directly applicable only to other FAST spectra. For the FAST velocities, the 
instrumental velocity offset with respect to the night sky lines is less than 10 km s~^ and is 
indistinguishable from zero (Kurtz Sz Mink 1998). We have found that the FAST template 
spectra are reliable also for the MMT and Hydra spectra; for those objects observed with several 
spectrographs, the redshifts derived from the separate spectra all agree within the errors. 

Many of the Hydra spectra have very low signal-to-noise ratios. We assess the robustness of 
the redshifts as follows: for each object, we construct all 2^ — 1 possible (non-zero) combinations 
of spectra from the N individual exposures. A'^ is 4 or 5 for most objects, and 8, 9 or 10 for 
objects observed in two configurations. For the few objects observed in 3 configurations, with 
N as large as 13, we examine only a subset of the 8191 possible combinations. We measure the 
redshift independently from each of these combined spectra. We accept a redshift for a target 
only if the following criteria all hold: (1) The i?-value increases with the number of individual 
spectra contained in the combined spectrum; (2) the velocity of the combined spectra stabilizes as 
the number of constituent individual spectra increases; (3) the best combined spectrum has an 
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i?- value >3. 

The comparatively narrow spectral range of the Hydra observations means that Ha is shifted 
off the edge of the detector for cz ^ 20000 km s^^. Cross-correlating emission line galaxies at 
these redshifts with the standard emission-line template results in a spurious redshift because 
the strongest remaining line in the spectrum is usually fit to Ha. We solved this problem by 
constructing a new template that differs from the standard emission-line template only by the 
absence of the Ha-nitrogen complex and the AA 6717,6731 A sulfur lines. 

We add an extra term of 65 km s~^ in quadrature to the errors of velocities determined 
from emission-line templates because the line emitting region may be offset in velocity from the 
systemic velocity of the galaxy (Kurtz & Mink 1998). 

Emission-line galaxies preferentially yield redshifts because an emission-line galaxy can be 
robustly matched to an emission template even when the signal-to-noise in the continuum is 
too low for a believable correlation with an absorption template. Thus, the completeness limit 
for emission-line galaxies is effectively fainter than for galaxies with no emission. The Hydra 
configuration observed through cirrus yielded few redshifts, resulting in a brighter effective 
magnitude limit; this effect is minor because there is considerable overlap on the sky among the 
four configurations. 

3.1. Spectral Classification 

We divide our spectra into two broad classes, emission (Em) and absorption (Ab), based 
on the presence or absence of Ha in the spectrum. Ha is "present" if its equivalent width 
exceeds 5 A. The robustness of this criterion is clearly a function of the signal-to-noise ratio of the 
spectrum; galaxies with low signal or with intrinsically weak Ha emission may be misclassified 
as absorption-line galaxies if the Ha peak does not rise appreciably above the noise. The 
misclassification is all in one direction: although some emission line galaxies may be misclassified 
as absorption, no absorption-line galaxy will be mistakenly classified as having emission. 

4. The Data 

Figure || (upper panel) shows the velocity distribution of all the observed galaxies. The 
cluster AWM 7, defined as in Koranyi et al. (1998) to consist of galaxies with redshifts in the 
range 2500-7500 km s^^, is cleanly and unambiguously distinguished from the background. 
The lower panel shows the velocity histogram of the cluster galaxies only; the solid histogram 
indicates the distribution of the entire sample, and the shaded histogram indicates the Em 
galaxies alone. Our cluster sample consists of 179 galaxies; we classify 42 as Em and 137 as Ab. 
A Kolmogorov-Smirnov test indicates that the probability of the Em and Ab samples being so 
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discrepant if drawn from the same distribution is small (PoyDahs ~ 0-0004); we conclude that the 
underlying distributions are different at the ~4(T confidence level. 

We compute the mean velocity and dispersion for the different galaxy populations following 
the prescription of Danese, De Zotti, & di Tullio (1980); the results are in Table |^. The Em 
galaxies have a much larger dispersion (~1100 km s~^ vs. ~600 km s^^), and are systematically 
offset in redshift by ~400 km s~^ from the Ab population (4800 km s~^ vs. 5200 km s~^). 
The extrema of the Em galaxies' velocity histogram are almost equidistant from the Ab mean 
velocity, but near 3800 km s-Hhere is an excess of Em galaxies above the otherwise almost 
uniform distribution; this excess decreases the mean redshift. Neither the 17 Em galaxies with 
3500 < cz < 4500 km s~^nor the 12 with 3700 < cz < 4300 km s~^ evince any strong clustering on 
the sky. Absent independent distance measurements to these galaxies, it is difficult to determine 
whether they are foreground interlopers, or whether the cluster Em galaxy population is indeed 
offset in velocity. 

The cD galaxy near the cluster center, NGC 1129, has a redshift of 5288 it 71 km s~^, which, 
within the errors, is at rest with respect to the absorption-line galaxy population. The peak of 
the ROSAT PSPC hard X-ray emission map, defined as the center of the pixel containing the 
greatest flux, is 6" from the nominal position of NGC 1129. We therefore adopt its position, 
{a,6)j2000 = (02:54:27.50, -^41:34:42. 50), as the center of the cluster. 

Figure |2| shows the distribution of galaxies with measured redshifts on the sky. Open 
circles and filled triangles denote Ab and Em galaxies in the cluster, respectively; crosses denote 
foreground or background galaxies. We do not plot targets for which we were unable to measure 
a redshift. The enhancement in the surface number density near the center of the cluster is due 
to the fainter magnitude limit reached by the MMT and Hydra observations; most Hydra targets 
turned out to be background sources. 

The Em galaxies are distributed much more uniformly on the sky than the Ab galaxies; 
the projected surface density of the Ab galaxies increases strongly toward the center, but the 
Em galaxies show no such enhancement. A two-dimensional Kolmogorov-Smirnov test (Press et 
al. 1992) indicates that the probability of drawing two such discrepant samples from the same 
underlying distribution is less than 5 x 10~^. The distribution of Ab galaxies is strongly flattened; 
there is a pronounced elongation in the east-west direction. 

The X-ray emission from AWM 7 is similarly elongated; Mohr et al. (1995) compute an 
axial ratio of 0.665 it 0.089 at a position angle of —83° it 7° from the X-ray flux within 7'.5 of 
the cluster center, in perfect agreement with the value of 0.67 computed by Dell'Antonio et al. 
(1995) for an isothermal ellipsoid model density distribution. Mohr et al. (1995) report X-ray 
ellipticities for 58 clusters; AWM 7 is more flattened than both the mean and median clusters, 
but there are 10 clusters that are more flattened still. In our case it is difficult to measure the 
flattening quantitatively from our optical sample for comparison with the X-ray results because 
the completeness limit varies across our survey. 
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Redshifts of cluster galaxies plotted against distance from the cluster center define a 
trumpet-shaped locus delimited by caustics (Kaiser 1987). Galaxies inside the caustics are bound 
to the cluster; those outside are unbound. We construct such a plot in Figure ^ (top panel); again, 
squares are Ab galaxies and triangles are Em galaxies. Although one can discern by eye where the 
caustics probably lie, they are in practice poorly defined. Their accurate determination requires 
many more galaxies than our sample contains (Diaferio 1999, Diaferio & Geller 1997, Geller, 
Diaferio & Kurtz 1999). 

From the velocity histograms and the position-velocity diagrams it is clear that the Em and 
Ab galaxies are not equivalent tracers of the cluster mass. A friends-of- friends algorithm (Huchra 
& Geller 1982) for finding substructure suggests that many of the Em galaxies are contained in 
smaller subgroups distinct from the main cluster (not shown); these may be accreting onto the 
cluster now, and, not yet being virialized, do not follow the dark matter distribution as closely 
as the Ab galaxies. There is considerable radial segregation of Ab and Em galaxies; the radial 
distribution of Em galaxies is more uniform and has a greater mean and median distance from 
the cluster center than that of the Ab galaxies. This behavior is observed in other clusters 
(Mohr et al. 1996, Adami et al. 1998), in groups (Mahdavi et al. 1999), and in the field as 
the morphology-density relation (Dressier 1980). Moreover, the velocity dispersion of the Em 
galaxies remains large with increasing radius while that of the Ab galaxies declines (see §5). These 
considerations suggest that the Em galaxies are not virialized tracers of the underlying dark 
matter distribution. Because our goal is an accurate determination of the mass profile, we exclude 
the Em galaxies from further analysis, following Mohr et al. (1996). We will revisit the issue of 
substructure in and around AWM 7 and other poor clusters in a subsequent paper. 

Figure ^ (middle panel) shows the velocity dispersion profile constructed only from the Ab 
galaxies. We sort the galaxies by distance from the cluster center, then evaluate the velocity 
dispersion in a sliding bin containing 11 galaxies; we plot the dispersion as a function of the mean 
distance to the galaxies in the bin. Neighboring points are strongly correlated, as they share 10 of 
their 11 galaxies; every 11th point (distinguished by error bars) is independent. 

Two large excursions from the overall trend dominate the profile. Two outlier galaxies, one 
at projected distance R ~ 0.7 Mpc and low velocity, the other at i? ~ 1.4 Mpc and high 
velocity (circled in the top panel of the figure) produce these excursions. It is clear that these two 
galaxies are responsible because the number of data points in each excursion equals the width of 
the sliding bin used to compute the profile. We therefore define a subset of the Ab sample, Abiss, 
excluding these two galaxies. We shall henceforth refer to the entire sample of 137 Ab galaxies as 
Abi37. The velocity dispersion profile for Abiss appears in the lower panel; the overall dispersion 
of these 135 galaxies is 543I3]' km s~^within R = 1.7 Mpc. After some fluctuations near the 
core, the profile becomes smooth and shows the characteristic declining profile seen in poor groups 
by Mahdavi et al. (1999) and in some A^-body simulations (e.g. Crone, Evrard & Richstone 1994). 
The decline in dispersion at small radii matches that of the "type C" clusters of Girardi et al. 
(1998). We revisit the issue of the cold core in §6. 
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5. Mass Profile 

We compute the mass profile of AWM 7 in two largely independent ways. First, we assume 
that the mass profile follows the functional form predicted by the hierarchical clustering simulations 
of Navarro, Prenk k, White (1995,1996,1997). This model has two parameters which we derive 
separately from the observations under the assumption that the Ab galaxies are virialized and that 
they trace the dark matter distribution. While this is a restrictive assumption, there is growing 
evidence that the NFW profile is applicable beyond the hierarchical clustering regime in which 
it was proposed. Huss, Jain &: Steinmetz (1999) find from high-resolution A^-body simulations 
that virialized dark matter halos with different formation histories nevertheless evince similar 
density and velocity dispersion profiles, and that this "universal" density profile is well described 
by the functional form of the NFW model. The collapse history of the cluster seems immaterial 
in determining the ultimate cluster profile. More concretely, Gcllcr, Diaferio &: Kurtz (1999) show 
that this profile describes the Coma cluster well for i? ^ 10 Mpc. For comparison, we also fit 
a non-singular isothermal sphere to AWM 7. 

We also compute the mass profile from the virial theorem alone. At small radii, this method 
gives incorrect results if not corrected for the exclusion of the outer cluster regions from the 
computation. The mass profile can be corrected with a surface term, which depends on both the 
mass and velocity dispersion profiles. These last two profiles are linked to the velocity anisotropy 
profile through the Jeans equation (e.g. Binney & Tremaine 1987), but the data are insufficient 
to constrain the problem completely. Usually one makes some assumptions about one or more of 
the profiles and solves for the others (e.g. Binney & Mamon 1982, The &: White 1986, Binney & 
Tremaine 1987, Girardi et al. 1998). Here, we make these assumptions only for the computation 
of the surface term correction to the virial mass estimator. 



5.1. The NFW model 



5.1.1. Description 
The spherically symmetric NFW density distribution is 

p(r) 

resulting in the enclosed mass profile 



^cPc 



{r/rc){l + r/rcf 



M{<r) = AixdcPcT, 



log(l -I- r/vc 



l + r/vc 



(1) 



(2) 



where Pc is the critical density of the universe, and 5c and are the two parameters of the model, 
corresponding to the overall normalization and the core radius, respectively. The radial scale in 
the context of cluster simulations is often quoted in units of r200) the radius where the mean 
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cluster density has dropped to 200 pc- NFW claim, on the basis of numerical experiments, that 
this radius approximately separates the virialized and infall regions (see also Cole &: Lacey 1996). 
NFW define the "concentration parameter" c = r2oo/^c- A low value of c arises if the core radius is 
comparable with the overall extent of the cluster, that is, the mass profile is not very concentrated 
toward the center. NFW demonstrate that more massive clusters are less concentrated, with the 
concentration parameter related to the overall normalization by 

^ 200 c3 
' 3 log(l + c) - c/(l + c) ■ ^ ^ 

Because the enclosed mass of the NFW profile diverges logarithmically, it is customary to quote 
M200 = 200/9c(47r/3)r2oo as the cluster mass. 



5.1.2. Fitting the NFW Model 

We can derive rc observationally from the distribution of projected radii of cluster galaxies. 
The projected surface number density of galaxies in the NFW model is given by 

ex l^li^ (4) 

where R = R/r^ is the projected radius in units of the core radius, and 

Yfj^. _ j {1- li^y^/^sech-^R ifR<l 
^^^^-\(R^-l)-y^sec-^R ifi?>l. 

Given the distribution of projected radii in our sample, we fit Vc with a maximum-likelihood 
technique by finding the value of that maximizes the probability 



Ri 



(6) 



of observing that particular distribution of projected radii, where Wi is a weight assigned to each 
galaxy to account for the non-uniform magnitude limit of our survey. We assign each galaxy one 
of two possible weights as follows: galaxies within the central 0.51 h^^ Mpc of the survey, where 
it is 97% complete to mn = 17.2, have w = 1.00; galaxies with Ri > 0.51 h^^ Mpc, where the 
survey is complete to m/j = 15.5, have w = 1.95, where 1.95 is the ratio of the number of galaxies 
with mji < 17.2 to the number with niji < 15.5 within 0.51 Mpc of the cluster center. This 
correction assumes that the luminosity function is the same in the core as in the outskirts of 
the cluster. We exclude galaxies fainter than mn = 17.2 within 0.51 Mpc and fainter than 
mji = 15.5 otherwise from the fit; there remain 53 galaxies with R < 0.51 h^^ Mpc and 70 
with R > 0.51 Mpc. Both outliers remain in the sample. The denominator in Eq. n serves 
to normalize the intrinsic probability of observing a particular galaxy i at projected radius Ri 
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(expressed in the numerator) to the possible range of R where it could be included in the sample. 
The upper limit of the integral is a function of rc, Rma.x{i'c) = Rwm/^c-, where Rw-ca is the limiting 
projected radius of the survey. In practice, we use the largest projected radius in the sample for 
i?iini- Table ^ lists the best-fit values of Tc for the Abias and Abi37 subsamples, computed from 
121 and 123 galaxies, respectively, along with 1-, 2-, and 3-c7 uncertainties, as determined by 
A(log£) = 0.5, 1.0, 1.5. There is little difference in the best-fit Vc for the two samples; if the Em 
galaxies are included, the best-fit Vc increases by ~50%, because the Em galaxies are less centrally 
concentrated (Fig. |3|). We compare the NFW model to the non-singular isothermal sphere in 



Section 5.2 



We use Monte Carlo simulations to assess the robustness of this method: for selected values 
of Tc, we construct 2000 sets of 123 projected radii drawn from the distribution of projected radii 
for an NFW profile with that Vc, we reproduce the non-uniform magnitude limit of our survey 
by only accepting galaxies with R > 0.51 Mpc with probability 1/1.95. We compute the 
maximum-likelihood recovered rc using Eq. ^ Figure Q shows the results for Tc in the range 0.04 - 
0.30 Mpc. The median recovered rc is within 0.002 of the input rc in all but one case; however, 
the distributions of the recovered rc are broad. The distribution of recovered rc is narrowest for 
small Tc, and broadens as rc increases, because (1) the peak in the probability distribution of 
projected radii becomes less pronounced as rc increases, and because (2) we sample a smaller 
dynamic range of R as rc increases (because the limiting projected radius Rum is constant). The 
distribution of the ratio of recovered to input rc is largely independent of input rc, however. The 
lower and upper la bounds (defined as the values of rc enclosing 34.2% of the distribution on 
either side of the median) are typically displaced —25% and -|-30% from the true rc, respectively. 
Increasing the sample size in the simulations to 500 radii (not shown) reduces the spread in the 
distribution of recovered rc by a factor of ~2. 

Knowing rc, the NFW model would be completely specified if we also knew r2oo, but this 
quantity is difficult to derive observationally. We can arrive at an estimate if we assume that the 
virial theorem holds, so that My = 3o"pr^G~^, where G is the gravitational constant, r„ is the 
virial radius, and dp is the global line-of-sight cluster velocity dispersion. Then 

- 100^0 (7) 

where we have used Pc{z) = 3Hq{1 + 2;)^/(87rG) for ^7 = 1, and the definition of M2oo- The 
right-hand side of this equation is observationally determined, and substituting for and M200 
using Eq. the left-hand side becomes a function of r^, r^, and r2oo which we solve for r2oo- We 
fit rc from the distribution of projected radii, and approximate ry by 

^ A^(iV-l) 

ry — — (8) 

where is the number of galaxies in the system, and Rij is the projected separation between 
galaxies i and j (Binney & Tremaine 1987). This estimator yields 1.30 h^^ Mpc for the 
Abi37 sample, and r^ ~ 1.28 Mpc for the Abiss sample. 
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The value of r2oo that solves Eq. ^ is fairly insensitive to the value of r^, in the regime of 
interest, where 500 < dp < 600 km s~^, r2oo varies ~3% as Tc ranges from 0.12 to 0.30 Mpc 
for any constant dp. The derived r2oo is more sensitive to dp, varying ~15% over the range 
500 < Up < 600 km s~^with constant Vc- For the obtained from Eq. ^, the values of r2oo are 1.02 
and 0.97 Mpc for the Abisy and Abias samples, respectively (Table |^. The corresponding 
value of the concentration parameter c, 4.6 in both cases, falls in the middle of the range 2.3-7.7 
seen by Carlberg et al. (1997) for rich systems, although the errors are large in both cases. 
Comfortingly, the values are also smaller than the c ~ 7.5 that Mahdavi et al. (1999) find for 
their lower-mass poor groups. For AWM 7, M200 is 2.47 and 2.13 x 10^^ h~^MQ for the Abisy and 
Abi35 samples, respectively. Here the larger mass is associated with a formally slightly smaller 
value of c (4.56 vs. 4.61), in keeping with the general trend seen by NEW, but this consistency 
may be accidental, considering the uncertainties in the computed values of c. It is interesting to 
note that two outliers (out of a sample of 137) can alter the recovered mass by ~15%. 

Because the cluster extends beyond r200 5 M200 in fact underestimates the total mass. 
Extrapolating the NEW mass profile to 1.7 Mpc (the maximum radial extent of the 
Ab sample) with Eq. || yields 3.49 and 3.01 x 10^^ /i^^M© for the Abisy and Abiss samples, 
respectively. 



5.2. Comparison to the Non-Singular Isothermal Sphere 

The non-singular isothermal sphere (NSIS) arises from setting finite boundary conditions at 
zero radius to the same differential equation that gives rise to the singular isothermal sphere (e.g. 
Binney &: Tremaine 1987, §4.4). It is characterized by a density normalization po and a core radius 
ro at which the the projected density falls to 0.5013 of its central value. In Eig. |5| we compare 
the fit of the projected surface number density from the NSIS and NEW models to the data. 



Although the maximum- likelihood technique of §5.1.2 makes optimal use of the data, it provides 
no independent assessment of goodness-of-fit. Here, although with some loss of information, we 
bin the observed surface number density for comparison with the models; we retain the same 
weighting scheme as in the maximum-likelihood calculation. 

The left side of Eig. ^ shows the observed profile (points with Poisson error bars), the NEW 
fit (solid line), and the NSIS fit (dotted line) for both the Abiss (top panel) and Abias (bottom) 
samples. The right side of the figure illustrates how the best-fit values of Tc, tq, and x^/^ vary with 
the choice of binning; the profiles and fits are shown for 11 galaxy equivalents per bin (galaxies 
beyond 0.51 Mpc count as 1.95 equivalents). 

Eor different binnings the quality of the fits of both the NEW and NSIS models, as measured 
by x^/^) varies much less for the Abias sample than for the Abi37 sample. The best-fit values of Vc 
and ro are quite stable for both samples. The decreased sensitivity to the binning provides further 
justification for excluding the two outlier galaxies. Eor both samples, the NEW model fits with a 
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lower jv than the NSIS model, independent of the binning (except in one case). The superiority 
of the NFW model is more consistent for the Abias sample, although this superiority is not 
overwhelming; for Abi35, NFW fits with 'x^ jv = 1.075, whereas the NSIS fits with 'x^ jv = 1.233. 
For Abi37, the respective x^/^ for NFW and NSIS are 1.184 and 1.252 (all for 11 galaxies per 
bin). There is thus some formal statistical justification for using the NFW model over the NSIS, 
but an unambiguous determination would require a survey out to 2 or 3 Mpc (cf. Geller, 
Diaferio & Kurtz 1999). 

With 11 galaxies per bin, the best-fit values of Vc are within 5% of the maximum- likelihood 
values for both samples, and well within the 1-cr range for all binnings, indicating that the loss of 
information inherent in the binning is not severe. 



5.3. The Virial Mass Profile 

We now compute the mass profile of AWM 7 from the virial estimator appropriate for the 
case of galaxies embedded in a diffuse distribution of dark matter, under the added assumption 
that light traces mass. The appropriate estimator (Binney &: Tremaine 1987) is 



2G E,;Ej<i|Ri-Rj|-i GN-1 ^' 



where Vi is the radial velocity relative to the cluster mean, and Rj is the projected radius vector 
relative to the center of the cluster. This estimator assumes that the galaxies are in dynamical 
equilibrium within the cluster potential, that galaxies trace the dark matter distribution, and 
that the entire cluster has been observed. These assumptions do not apply to the Em galaxy 
population; indeed, when we apply this prescription to the Em galaxies alone, the resultant mass 
estimate is a factor of ~4 greater than for the Ab galaxies. To correct for the failure to include 
the entire angular extent of the cluster in our sample, we correct the mass estimate with a surface 
term C{b) (The &; White 1986) that depends on the limiting radius b of the observations. The 
corrected mass is 



M,^(< h) = Mest[l - C(6)] = Mest 



/q 47rr^p(r) dr C7^(< b) 



(10) 



where a{< b) denotes the integrated velocity dispersion within the limiting radius 6 (e.g. Girardi 
et al. 1998). 

The surface term incorporates information about the velocity dispersion profile and the 
density profile; the former is measurable from the data, but some assumption must be made about 
the form of the density profile — in order to correct that very profile. The &: White (1986) and 
Girardi et al. (1998) posit a functional form for the mass profile, which they then fit to the data 
on the basis of the projected number density of galaxies. Here we assume an underlying NFW 
mass profile. 
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5.3.1. Behavior of The Surface Term 
The radial dependence in the NFW model of the surface term in Eq. [in is 



Cnfw(?') 



u 



2 r -1-1 r _ -i2 



[l+U) 



log(l + u) 



u 


-1 




l + u_ 




.o-(< r)_ 



(11) 



where u = r/rc. The last term involves both the overall velocity dispersion a and the radial 
velocity dispersion ar , and is thus a function of the anisotropy parameter /? = 1 — (t|/(T^ , which is 
itself a function of r. Now o"^ = o"^ + (t| + o"^; under the assumption of spherical symmetry the 
last two terms are equal, so a = {3-2p)a^, which we substitute into the denominator of the last 
term for computation. The anisotropy and velocity dispersion profiles are not independent; they 
are coupled via the Jeans equation. Girardi et al. (1998) find that for their mass model a constant 
anisotropy profile /3(r) = produces a velocity dispersion profile that peaks at ~0.1 r„, dips 
slightly in the core, and decreases gradually at large radii; the simple functional form /3(r) = —k/r, 
where A; is a constant, produces a velocity dispersion profile with a more severe dip in the core. In 
this case, the orbits in the core become purely tangential as /3 — > —oo. 

Figure ^ shows the behavior of the surface term for an NFW mass distribution and the 
observed velocity dispersion profile (Fig. ^) under the assumption of three different anisotropy 
profiles: constant (top panel), P{r) = —k/r (middle), and /3(r) = 1 — k/r (bottom). The multiple 
curves in each panel correspond to the values of /3 (top) or k (middle and bottom) enclosed in 
brackets. One of the extreme curves is labeled in each case for reference. The values of k are 
chosen to bracket the value inferred from the profiles in Fig. 3 of Girardi et al. (1998). Note that 
the three panels are not drawn to the same scale. Because these classes of /? profile are derived 
from fits to the mass model of Girardi et al. (1998), they are not a priori strictly valid for the 
NFW profile; nevertheless, it is reasonable to adopt them here as representative of the general 
forms that the profiles may have. 

By definition, the surface term must be everywhere less than unity (to ensure a positive 
mass). The results shown for constant large values of /3 are thus clearly unphysical. At small radii, 
the data are too noisy for all our assumptions to hold; a proper treatment would require a full 
solution of the Jeans equation. It is unlikely that this treatment would improve the results very 
much, however, because of the intrinsic noisiness and discreteness of the data. At larger radii, the 
variations in the surface term become less pronounced, and the surface term in general becomes 
less important, decreasing to ^ 5% for each profile at the outer limit of our sample. Sufficiently 
far from the core, then, the virial mass estimator provides a robust estimation of the enclosed 
cluster mass, but at small radii the surface correction is too variable to result in a reliable mass 
profile. The virial mass estimator thus seems to be a good tool for estimating total cluster masses, 
but it is not well suited to the determination of mass profiles in cluster cores, where the surface 
term correction is most important but least well determined. 



A reliable description of the surface term at small radii may only be possible from cluster 
simulations in which the dark matter density and velocity distributions are fully known to the 
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experimenter. In that case, empirical surface term profiles can be built from the data and 
compared to reconstructions that an observer would generate under various assumptions of the 
anisotropy profile. 

5.3.2. The Virial Mass Profile 

Figure ^ shows the virial mass profile (Eq. with no surface term correction, as derived 
from three subsets of our galaxy sample: the entire sample (crosses), the Abi37 sample (triangles), 
and the Abiss sample (squares). The entire (Em+Ab) sample is included only to illustrate that 
the inclusion of the emission-line galaxies with their larger velocity dispersion skews the results 
towards higher mass. Indeed, this difference is evident in the two Ab profiles as well: they agree 
up to the first outlier, at which point the Abisy mass profile jumps; the two Ab profiles then grow 
more or less in parallel until the next outlier. We compute uncorrected enclosed masses within 1.7 

Mpc of 3.10 and 2.74 x 10^^ h~^MQ from the Abi37 and Abias samples, respectively. The 
surface term correction at this radius is on the order of 2-5%, depending on the anisotropy. 

We estimate the error in the enclosed mass profile by the statistical bootstrap method 
(Diaconis &: Efron 1983): for the N galaxies enclosed within a given projected radius, we compute 
the possible enclosed masses from all subsets of — 1 galaxies. We take the standard deviation 
of the resultant set of masses as the formal error in the estimated enclosed mass. The formal 
errors become quite small, on the order of 3%, at large radii. When the number of galaxies is large 
and the distribution of velocities is well-behaved, the random rearrangement of velocities has little 
effect on the mass estimate; thus, the distribution of the bootstrapped masses is narrow. The 
actual error is in fact larger, as demonstrated by the sensitivity to the outliers. The error in the 
enclosed virial mass is ~15%. 

The solid smooth curves indicate the NEW mass profiles we fit previously, and the dotted 
lines indicate the profiles corresponding to the la ranges of Vc from Table ^ with r2oo recomputed 
from Eq. |^ The agreement with the virial mass is quite good, to within ~15% at 1.7 Mpc for 
both Ab samples, with the virial mass the smaller of the two in each case. The agreement along 
the profile in general is reasonable given the uncertainty in the surface term. The discrepancy is 
greatest near R ~ 0.5 Mpc, corresponding to r/rc ~ 2, where the surface term correction 
can be on the order of 20% in the proper direction; at the virial radius the masses are in perfect 
agreement, and at larger radii the NEW mass profiles overtake the virial mass profiles. Because 
r2oo is insensitive to Tc, the variation in Tc results in little difference in the mass profile at large 
radii, but of course has a greater effect on scales of a few times r^. Table |5| compares the NEW 
and virial masses. Based on the sensitivity of the virial mass to the outliers and on the effect of 
the surface term, we assign an error of 0.5 x 10^^ Mq to our mass estimate. 

Our mass profile is in agreement with the ROSAT X-ray mass of Dell'Antonio et al. (1995), 
who derive M(< 0.27) = 0.8 x 10^^ M© (indicated by the arrows in Eig. |^), and is within the 
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allowable range of mass profiles determined from ROSAT data by Neumann &; Bohringer (1995) 
and from ASCA data by Loewenstein & Mushotzky (1996) and Markevitch & Vikhlinin (1997). 
The X-ray masses arc all quite uncertain and extend only to i? 0.5 h"^ Mpc. Koranyi et al. 
(1998) illustrate the range of mass profiles from X-ray determinations; the range is so broad as to 
include even the virial mass profile computed from emission-line galaxies; any plausible kinematic 
mass profile is thus likely to be consistent with the current X-ray limits on the mass. The increased 
sensitivity of the Chandra Observatory will allow X-ray mass profiles to be determined with more 
precision and to larger radii. 

6. Discussion 

Galaxy clusters display a range of velocity dispersion profiles ( e.g. den Hartog & Katgert 
1996, Girardi et al. 1998). Some are peaked at the center, others are more or less flat, still others 
rise with increasing radius near the center. This variety suggests that clusters today are still in 
a variety of dynamical states, and that indiscriminately combining a set of clusters to improve 
the statistics can easily lead to incorrect results if clusters with different kinematics are averaged 
together. Even within a single cluster, it is important to discriminate between the generally red, 
virialized population and the newly infalling bluer galaxies; including the latter in steady-state 
kinematical analyses will artificially inflate the mass estimate; we suggest that this effect accounts 
for the mass of 5.77t\io h'^M^ derived for AWM 7 by Girardi et al. (1998). 

In principle, the non-virialized Em galaxies can be used as tracers of the escape velocity to 
arrive at an independent mass measurement (Diaferio &: Geller 1997, Diaferio 1999). In practice, 
there are not enough Em galaxies in our sample for this method to yield a robust mass estimate; 
although the formal mass returned by this method for AWM 7 (kindly computed for us by A. 
Diaferio) is comparable to the ones we derive here, the fractional error in that mass is ~0.9. 

Our enlarged sample confirms the presence of a cold core in the cluster (Koranyi et al. 1998); 
although still noisy near the cluster center, the velocity dispersion within a projected radius of 
~0.2 Mpc is lower than in the range 0.2-0.5 Mpc. Beyond ~1 h^^ Mpc the velocity 
dispersion drops below the range in the core. Abell 576 (Mohr et al. 1996), a similarly well-sampled 
cluster, has a more pronounced cold core but a rising absorption-line velocity dispersion profile in 
the range 0.4-1.2 Mpc. A simple dynamical argument suggests that cold cores arise naturally 
from the NFW density profile. Assuming virial equilibrium to hold in the core, cr^ oc GM{< r)/r. 
Modeling the density profile at small radii as a power law p{r) oc r~" yields M(< r) oc r^~", so 
a oc 

^i-a/2 p^j. ^j^g -^YW profile a = 1 in the core, resulting in a rising velocity dispersion profile. 

Because the velocity dispersion profile, mass profile, and orbital anisotropy profile are 
coupled through the Jeans equation, it is difficult to derive a mass profile without making some 
assumptions about the form of one or more of these functions along the way. The projected 
velocity dispersion profile is directly observable, but it is sensitive to outliers and substructure, and 
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requires either extensive smoothing or a large number of galaxies for a robust determination. The 
anisotropy profile is also difficult to constrain, with some evidence for a variety of profiles possible 
in clusters (Girardi et al. 1998). The mass profile can be estimated from X-ray observations, 
but these tend to have few data points, large errors, and cover only the inner regions of clusters; 
consequently, they can be fit by a wide range of models, resulting in large uncertainties in mass 
when extrapolated to large cluster radii. 

It seems difficult to avoid circular logic in computing the surface term, as it depends on both 
the anisotropy profile and the density distribution that one is trying to determine. One can make 
a self-consistency argument, as we do here, to verify that the assumed density profile used to 
correct the virial mass profile agrees with the virial profile so corrected; if the surface term were 
small, this circularity would be less troublesome, because the assumed density profile would affect 
the virial mass profile only through that small correction. But as Fig. ^ illustrates, the surface 
term correction is often significant, particularly at small radii. Not only is the correction large, but 
it is also uncertain, sensitive to uncertainty in the anisotropy profile alone. The coupling between 
the mass profile and the surface term is thus quite strong; even with the assumption of a density 
profile (which in principle renders the whole exercise redundant because the mass profile is then 
already known), the surface term correction is too poorly constrained to produce a reliable mass 
profile. Only at the edges of the cluster does the surface term correction become reliable, but 
by then it is no longer very important. Accurate determination of cluster masses from the virial 
estimator, then, requires broad angular coverage; dense sampling of the core is not enough. 

7. Conclusion 

We derive the mass profile of AWM 7 in two independent ways from an optical sample 
of 179 galaxies with redshifts within 1.7 h^^ Mpc of the cluster center; we use only the 137 
absorption-line galaxies (and a 135-element subset of them) in the analysis. One method uses 
the surface number density and global velocity dispersion; the other uses individual velocities and 
pairwise distances. Both assume and underlying NFW profile. The two methods yield remarkably 
consistent mass estimates, particularly at the outer limit of our sample; the enclosed mass within 
a projected radius of 1.7 h^^ Mpc agrees within 15% for both samples. The mass estimates are 
sensitive to the inclusion of two outliers in the sample; both the NFW and virial masses increase 
by ~15% when they are included. Including the Em galaxies, with their larger dispersion, in the 
calculation results in a virial mass larger by a factor of two. The cluster mass within 1.7 Mpc 
is ~(3±0.5) X 10^"^ h-^MQ. 

The surface term correction to the virial mass estimator is sensitive to noise in the velocity 
dispersion profile and to the assumed velocity anisotropy profile. The surface term at small radii 
is too unreliable for the virial mass estimator to be used for measuring mass profiles; at large 
radii, where the surface term is more certain, it is no longer an important correction to the 
enclosed mass. The virial mass estimator is thus reliable for measuring total cluster masses, but 
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not enclosed mass profiles; to overcome the limitations of the surface term at small radii, cluster 
redshift samples must have broad angular coverage. Dense sampling of cluster cores can reduce 
the scatter in the surface term if there are enough galaxies in the core, but our experience with 
AWM 7 (22 cluster members out of 154 Hydra redshifts) suggests that telescope time is better 
spent probing the outskirts of the cluster. 

Untangling the relation between the anisotropy profile, velocity dispersion profile, and density 
profile in cluster cores may ultimately require high-resolution simulations that track both dark 
matter particles and, once they have formed, galaxies. 
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spectra reported herein, and Susan Tokarz for preliminary reductions. We thank Mike Kurtz and 
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Table 1. Redshifts of Target Galaxies 
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Table 1 — Continued 
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Note. — RA and Dec (J2000), and redshift and error (in km s ^) for galaxies with measured 
redshifts. 
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Table 2. Velocity Distribution 



Sample 


Ngal 


cz 
[km s~-^] 


[km s~-^] 


dp 
[km s~^] 


^200 

[km s~-^] 


Em+Ab 
Em 

Abi37 
Abi35 


179 
42 
137 
135 


5186.2 

4821.5 
5247.6 
5256.6 


54.3 
168.0 
49.0 
46.7 


727.0t^^| 
1088.8;i^|^ 

c;7o 7+38.3 
^' -J-' -32.0 

542.9li;6 


601.51^^ 

578.ili:i 



Note. — Velocity dispersions are computed from galaxies within 
1.72 Mpcof the cluster center for Ab samples, and within 1.80 
hr^ Mpcfor Em and Em+Ab samples. (T200 is the projected velocity 
dispersion within r2oo (§5.1). 



Table 3. NEW from Projected Radii 



Sample ''c / [ ^ ^ Mpc] Icr range 2(7 range 3o" range 



Abi37 0.223 0.174 - 0.285 0.157 - 0.316 0.145 - 0.342 

Abi35 0.211 0.164 - 0.270 0.148 - 0.299 0.137 - 0.324 



Note. — Tc is computed from two-zone magnitude-limited samples of 123 and 121 members, 
respectively. 
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Table 4. NFW Profile Parameters 



Sample 






c 


5c 




M200 


M(< 1.7 h-^ : 




[/i-i Mpc] 


[/i-i Mpc] 






[/i-i Mpc] 




[IQi^/i-iM, 


Abi37 


0.223 


1.02 


4.6 


7.1 X 10^ 


1.30 


2.47 


3.49 


Abi35 


0.211 


0.97 


4.6 


7.2 X 10^ 


1.28 


2.13 


3.01 



Note. — Tc is computed from two-zone magnitude-limited samples of 123 and 121 members, 
respectively. 



Table 5. NFW and Virial Masses 



Sample 


Mnfw(< 1-7 h-^ Mpc) 


Mvt(< 1.7 /i"^ Mpc) 


AM/MvT 




[lO^^/i-iMo] 


[lO^^/i-^Mo] 





Abi37 
Abi35 



3.49 
3.01 



3.10 
2.74 



0.126 
0.099 
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Fig. 1. — Top: Velocity histogram of all observed targets. Bottom: Velocity histogram of cluster 
members. The solid histogram denotes the entire sample (both Em and Ab galaxies); the shaded 
histogram denotes Em galaxies only. The dotted vertical rule indicates the mean redshift of the Ab 
galaxies; the solid vertical dash over the histogram indicates the redshift (and error) of the central 
cD galaxy. 
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Fig. 2. — Position on the sky of all objects with measured redshifts. Note the denser sampling at 
the center, in the region covered by the Hydra and MMT observations. Crosses denote foreground 
and background objects. Open circles denote Ab galaxies in the cluster. Filled triangles denote 
Em galaxies in the cluster. Coordinate axes are in decimal degrees. 
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Fig. 3. — Top: Velocity as a function of projected distance from the cluster center. Ab galaxies 
are filled squares, Em galaxies are open triangles. The two outliers are circled. Middle: Velocity 
dispersion profile for the 137 Ab galaxies in the cluster sample. Every eleventh point is independent 
and is shown with la error bars. The horizontal rule indicates these galaxies' overall dispersion of 
574 km s~^. Bottom: Same as middle panel, but with two outliers removed. Horizontal rule is now 
at 543 km s~^. 
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Fig. 4. — Histograms of recovered values of the NFW parameter Tc from 2000 Monte Carlo 
realizations of sets of 123 projected radii drawn from NFW models with varying Tc- The histograms 
are all to the same scale, but for the highest input Tc may extend off the right edge of the plot. The 
triangles in the histograms and the numbers along the vertical axis both indicate the input value 
of re- 
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Fig. 5. — Comparison of the actual binned projected surface number density of galaxies to the 
NFW (solid) and non-singular isothermal sphere models (NSIS, dotted). The top half of the figure 
refers to the Abiss sample; the bottom half to the Abi37 sample. Left: the actual surface density 
profile is shown as the dots with error bars; there are 11 galaxies per bin. The curves show the 
best-fit NFW and NSIS models. Right: dependence of the fits on the binning. The behavior of the 
best-fit value of rc (for the NFW model) or ro (for NSIS) and of ~)(^ jv is shown as a function of the 
binning. The horizontal rules indicate the maximum-likelihood values of Vc for the two samples. 
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Fig. 6. — Behavior of the surface term correction to the virial mass estimator as given by Eq.|10 
assuming an underlying NFW mass profile for three models of the anisotropy profile P{r). The 
assumed NFW profile uses = 0.211 (the best-fit value for the Abias sample), and we use 
the observed velocity dispersion profile. Top panel: (3{r) = k; middle: /3(r) = —k/r; bottom: 
P{r) = 1 — k/r, where k is a constant in each case. There is one curve for each value of k listed in 
the label, with one of the extreme curves labeled for reference. 
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r/h~^Mpc 

Fig. 7. — Enclosed mass profile as determined from the virial theorem, for Em, Abi37, and 
Abi35 samples. Error bars show the formal bootstrap errors (see text). The solid curves are 
the independently derived NEW profiles tabulated in Table the dotted curves are the profiles 
corresponding to the la ranges of rc from Table |^ with r2oo recomputed from Eq. |^ The arrows 
point to the X-ray mass of Dell' Antonio et al. (1995). 



